Fermionic theory for quantum antiferromagnets with spin S > 1/2 
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The femiion representation for 5* = 1/2 spins is generalized to spins with arbitrary magnitudes. 
The symmetry properties of the representation is analyzed where we find that the particle-hole 
symmetry in the spinon Hilbert space of S =1/2 fermion representation is absent for 5* > 1/2. As 
a result, different path integral representations and mean field theories can be formulated for spin 
models. In particular, we construct a Lagrangian with restored particle-hole symmetry, and apply 
the corresponding mean field theory to one dimensional (ID) 5* = 1 and 5* = 3/2 antiferromagnetic 
Heisenberg models, with results that agree with Haldane's conjecture. For a S = 1 open chain, we 
show that Majorana fermion edge states exist in our mean field theory. The generalization to spins 
with arbitrary magnitude S is discussed. Our approach can be applied to higher dimensional spin 
systems. As an example, we study the geometrically frustrated S = 1 AFM on triangular lattice. 
Two spin liquids with different pairing symmetries are discussed: the gapped px + ipy-w&ve spin 
liquid and the gapless /-wave spin liquid. We compare our mean field result with the experiment 
on NiGa2S4, which remains disordered at low temperature and was proposed to be in a spin liquid 
state. Our fermionic mean field theory provide a framework to study 5* > 1/2 spin liquids with 
fermionic spinon excitations. 

PACS numbers: 75.10.Kt, 75.10.Jm, 71.10.Hf 



I. INTRODUCTION 

Quantum magnetism is one of the oldest central prob- 
lems in condensed matter and many body physics. Few 
exact results are known except at one dimension or when 
long range spin order is established at low temperature. 
In the later case. Landau's symmetry breaking paradigm 
is applicable, and the long range spin order determines 
the low energy physics of the system. The low energy 
spin excitations are spin waves or magnonsi System- 
atic corrections to spin wave theory can be obtained 
through semiclassical approach (1/5* expansion^) and/or 
nonlinear a model.'^ However, it was suggested that for 
some anti-ferromagnetic materials, strong quantum fluc- 
tuations due to small spin magnitude and low dimen- 
sionality combined with geometric frustration may lead 
to quantum coherent, spin disordered ground states.-*^'^ 
In such cases, a quantum fluid approach^ is probably a 
better starting point to describe the low energy physics 
instead of usual spin wave or semiclassical approach. 

There exist at present two common quantum fluid ap- 
proaches based on different "particle" representations of 
spins. In the "bosonic" approach^'^ spins are represented 
by spin up and down Schwinger bosons with the con- 
straint that total boson number at each site is 25. A 
spin-disordered state is obtained in a mean field theory 
as long as Bose condensation does not occur. The low ly- 
ing excitations coming from this theory are bosonic. The 
other possibility is to represent spins by fermions. Up 
to now, the fermionic approach is mainly restricted to 
S = 1/2 spin systems^'^-. Because of the Pauli exclusion 
principle, we cannot put more than two spin up or down 
fermions on a site to form 5 > 1/2 objects as in bosonic 
approach. One way out is to introduce more species of 
fermions. It was proposed in Ref. [Hi, that {2S + 1)- 



species of boson/fermions can be used to construct the 
spin-5 swapping operators. Later, it is shown that the 
fermionic representation can be used to describe a S' = 1 
valence bond solid state.— 

In this paper, we systematically study the fermionic 
representation for arbitrary spins. To construct a spin 
operator with spin magnitude S, we follow Ref. [ll| and 
introduce different fermion operators for different Sz 
states. We find that for half-odd-integer spins, the spin 
operator thus constructed is invariant under a SU{2) 
transformations^^ i.e. the symmetry group is SU{2) 
group whereas for integer spins the symmetry group is 
U{1)(S)Z2- Fermionic path integral formulations for spin 
systems with arbitrary spin-S* can then be formulated 
and the corresponding mean field theories can be stud- 
ied. We shall see that one of the mean field theories 
applied to ID antiferromagnetic Heisenberg model pro- 
duces mean field excitation spectrums which agrees with 
Haldane's conjecturei^ Moreover, there exist zero energy 
Majorana edge modes for open integer spin chains. We 
also apply our approach to two-dimensional (2D) anti- 
ferromagnetic Heisenberg models on triangular lattices 
where we obtain two translational and rotational invari- 
ant spin liquid solutions, one (with p^ -\- ipy-w&ve paring 
symmetry) is gapped and the other one (with /-wave 
pairing symmetry) has nodes at the spinon Fermi level. 
Comparing with experiments on NiGa2S4, we argue that 
the /-wave pairing spin liquid is a plausible ground state. 

It should be noted that rather different 
fcrmionic/bosonic spinon approaches have also been 
applied to SU{N) spin models where the dynamic 
symmetry is generalized from SU{2) to SU{N), and the 
Hamiltonian is composed of SU{N) generators i^'^d" . 
In these approaches, A^-species of particles are in- 
troduced to construct the SU{N)-spin operator and 



Hamiltonians.^- The major difference between our 
approach and the SU{N) approaches is that (25 + 1) 
species of fermions are introduced to form a S'C/(2)-spin 
(or irreducible representation of SU{2) group) in our 
approach. Unhke the SU{N) models^^ our approach is 
not hmited to bipartite lattices. 

The paper is organized as follows. In section II, we 
introduce the general fermionic spinon representation for 
arbitrary spin S and discuss the symmetry properties 
associated with the representation. In section III, we in- 
troduce the path integral formalism where different ways 
of imposing particle number constraints are presented. 
A particle-hole symmetric way of imposing constraint is 
discussed. In section IV. A and IV. B, we focus on ID 
Heisenberg models for spin-1 and spin-3/2, respectively, 
and study two versions of the mean field theories. The 
cases for general spin S are discussed in section IV. C. In 
section V we study 2D AFM models on triangular lattice. 
Our paper is summarized and concluded in section VI. 



II. GENERAL FERMIONIC SPIN 
REPRESENTATION 

We begin with the fermionic representation for spins. 
For spin S — 1/2, two species of fermions c^ and c^ are 
introduced to construct the three spin operators S^'^'^. 
To generalize this fermionic representation to arbitrary 
spin-S, we introduce 25* -1-1 species of fermionic operators 
Cm satisfying anti-commutation relations. 



\ C„j ,Cj^l — 0,, 



(1) 



where m, n — S, S — 1,- ■ ■ , —S. The spin operator can 
be expressed in terms of Cm and cj^'s, 

S = C^IC, 

where C = {cs,cs-i,- ■ ■ ,c-s)'^ and /° (a = x,y,z) is 
a (25* + 1) X (25* + 1) matrix whose matrix elements are 
given by 



ja 



{S, missis, n). 



It is easy to see that the operators 5*" satisfy the SU{2) 
angular momentum algebra, [S'°,5'''] = ie°'^'^S'^. Under a 
rotational operation, C is a spin-S "spinor" transform- 
ing as Cm — >■ D^^Cn and S is a vector transforming 
as S"- — >■ RabS'', here D^ is the 25 -I- 1-dimensional ir- 
reducible representation of SU{2) group generated by I 
and R is the adjoint representation. 

As in the S = 1/2 case, a constraint that there is only 
one fermion per site is needed to project the fermionic 
system into the proper Hilbert space representing spins, 
i.e. 



(7V,-iV/)|phy)=0, 



(2) 



where i is the site index and Nf = 1 (the particle rep- 
resentation, one fermion per site). Alternatively, it is 



straightforward to show that the constraint Nf — 2S 
(the hole representation, a single hole per site) repre- 
sents a spin equally. The Nf = 1 representation can be 
mapped to the Nf = 25* representation via a particle- 
hole transformation. For S — 1/2, the particle picture 
and the hole picture are identical, reflecting an intrinsic 
particle-hole symmetry of the underlying Hilbert space 
which is absent for S > I. 

Following Affleck et al. ^ we introduce another 

"spinor" C = (cls,-cls+i,cLs+2'--- .(-l)'H)^' 
whose components can be written as Cm — (~1)'^ ™cLm, 
where the index m runs from S* to —5 as in C. To 
examine whether C is really a "spinor" , we construct 
a spin singlet state for two spins at site i and j. 



V2S+1 I 



V2S+ 



hf^L-si-^r-ns,m},\s, 



Therefore CjCj is a scalar operator as CjCj, meaning 
that C and C must behave identically under spin rota- 
tion. Consequently, the spin operators can also be writ- 
ten in terms of C, 

S = C^IC. 

Combining C and C into a {2S + 1) x 2 matrix ifj = 
(C, C),® we can reexpress the spin operators as 



S = iTr(^ti^) 



(3) 



and the constraint can be expressed as 

TT{ija,^ip'')^2S+l-2Nf ^±{2S-l), (4) 

where the + sign implies Nf — 1 and — sign implies 
Nf = 2S.^ 

We are interested in two kinds of operations (or groups) 
acting on tp. The first kind belongs to a {2S + 1)- 
dimensional irreducible representation of SU(2) group, 
acting on the left of ip. Suppose G is an element of this 
irreducible representation, then G^ I°'G = R"^^!^ {R be- 
longs to the adjoint irreducible representation of SU{2) 
group, which is a iS'0(3) matrix). Under transformation 
-0 — >• Gip, the spin operator S"" — ^Tr {ip'^ I°-ip) becomes 

S"' — >■ R'^^S^, which means a rotation of the spin. It is 
obvious that the particle number constraint Eq. ^ re- 
mains unchanged under the action of G. 

The other kind belongs to a 2 x 2 unitary group acting 
on the right of ip which keeps the spin operator Eq. Q 
and the fermionic statistics Eq. ([T]) invariant. This 
group reflects the symmetry properties of the underlying 
Hilbert space structure in the fermionic representation. 
We call it an internal symmetry group. The internal sym- 
metry group is different for integer and half-integer spins. 
It is C/(1)(§)Z2 = {e^'^-^cr^e^'"'=^ = e-"'-^a^;e e R} for 
the former and SU{2) for the latter. We leave the rigor- 
ous proof in appendix A, and shall explain qualitatively 
the reason behind here. Notice that C and C are not 
independent. The operators in the internal symmetry 
group will "mix" the two fermion operators in the same 



-S+1' 



row of C and C, i.e. cs and c_g, cs-i and —c'_ 
etc. For integer spins, co and (— 1)"^Cq will be "mixed". 
To keep the relation {cq, cJ} = 1 invariant, there are only 
two methods of "mixing": one is an C/(l) transformation, 
the other is interchanging the two operators. These op- 
erations form the U{1)®Z2 group. For half-odd-integer 
spins, the pair (cq, (— 1)'^cJ) do not exist, and the symme- 
try group is the maximum group SU{2). The difference 
between integer and half-integer spins is a fundamental 
property of the fermionic representation as we shall see 
more later. 

Now let us see how the constraints Eq. ^ transform 
under the symmetry group. For S = 1/2, the constraint 
Eq. (J4]) is invariant under the transformation ip — >■ ^W 
because the right hand side vanishes (due to the particle- 
hole symmetry of the Hilbert space). For integer spins, if 
W = e*"^"^, then Wcr^W^ = a^-, and Eq. dU is invariant. 
\iW = a^e^'''^ then Wa:,W'^ ^ -o^, meaning that 
the "particle" picture(-f sign in Eq. ^) and the "hole" 
picture (— sign in Eq. ((4])) are transformed to each other. 

For a half-odd-integer spin with S > 3/2, W G SU{2) 
is a rotation and we need to extend the constraint into 
a vector form similar to 5 = 1/2 case^^ so that Eq. ((4]) 
becomes 



Tr(^a^T)^ (0,0, ±(25-1))^ 



(5) 



Under the group transformation tp — )■ ipW, 

Tr(V;aV;^) ^ [R-'M 0, ±{23 - l)f. (6) 

where Wa°'W'^ — Rabcr'', a,b — x,y,z, i.e. i? is a 3 
by 3 matrix representing a 3D rotation in the internal 
Hilbert space. The transformed constraint represents a 
new Hilbert subspace which is still a (2iV-|-l)-dimensional 
irreducible representation of the spin S'C/(2) algebra. Any 
measurable physical quantity such as the spin S remains 
unchanged in the new Hilbert space. Therefore, for 
half-odd-integer spins {S > 3/2), there exists infinitely 
many ways of imposing the constraint that gives rise to 
a Hilbert subspace representing a spin. However, for in- 
teger spins, there exists only two possible constraint rep- 
resentations. 

The different representations of constraints for S > 1/2 
systems result in different path integral representations 
and different mean field theories. These mean field the- 
ories are equivalent in the sense that they can be trans- 
formed to each other by the internal symmetry group. In 
the next section, we shall study the Heisenberg model in 
different representations and shall construct a "mixed" 
Path Integral representation which restored particle-hole 
symmetry, and the internal symmetry group becomes "al- 
most" a gauge symmetry.^ The "mixed" representation 
is studied in section IV where a new mean field theory 
is proposed which recovers Haldane conjecture at ID for 
the Heisenberg model. 



III. PATH INTEGRAL FORMALISM AND 

MEAN FIELD THEORY FOR HEISENBERG 

MODEL 

We shall focus on the antiferromagnetic Heisenberg 
model H = Jj^a j) ^i ' ^j with J > in the rest of 
the paper. We start by presenting some useful algebraic 
manipulations of the Hamiltonian, and then discuss the 
general path integral formalism. The different ways of 
handling constraints and how they affect the symmetry 
of the Lagrangian will be discussed in the process. 



A. Heisenberg model and an effective Hamiltonian 

It is known that for spin-1/2 the Heisenberg interaction 
can be written as^ 



a,; • a, 



-Tr : (V'lV'jV'] 

i • fv^ L At 

4 " 



{xlx^J + ^lj\j) ■; 



where 



X^J 



c}c,, 



cjc,. 



(7) 



(8) 



are spin-singlet operators. Xij ^^id Ay will be extended 
to arbitrary spin magnitudes with the above definition 
and will be used in the following discussions. Interest- 
ingly, an expression almost the same as Eq. ([7]) holds for 
5' = 1, but this is no longer true for larger spins. We 
shall present precise formula for S = 3/2 and 5 = 2, and 
provide a general discussion for higher spins. 
For S — 1, the three spin matrices are 



(I-V 



where I± = I^ ± ily The matrix operator ^p is given by 



V2 \ 




/l 


V2 


,/. = 





/ 




\o -1 



.p^icc) 



Ci c'_^ 



The spin operator can be written in form of Eq. ([3]) , and 
it can be shown after some straightforward algebra that 
the Hamiltonian can be written as 



H 



J 



J ^ S, . S, = -- ^ Tr : (VlV-.^jV'.) 
(id) (id) 

For S = 3/2, the three spin matrices are 

, = {I-y 



(9) 



fo Vs \ 




/f \ 


2 




i 


\/3 


,/. = 


-i 


\o y 




[o -I) 



and the matrix operator xp is given by -0 = 
iC C) where C = (c3 ,ci ,c_i ,c_3)^ and C = 

(c_3 , — c_j_, cl , — €3)"^. In this case the singlet operators 
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Xij and Ay alone are not enough to represent the Heisen- 

berg Hamiltonian, and triplet hoping and pairing terms 

are necessary. After some straightforward but tedious 

algebra, we find that the Hamiltonian for spin S = 3/2 

Heisenberg model can be written as 

H ^ J y Si • Sj 

(10) 



It is interesting to note that the above form holds also 
for 5* — 1^^ . However, for 5 — 1, the two terms in 
the square bracket can be transformed from one to an- 
other, which does not hold for 5* ~ 3/2. This suggests 
that S = 1/2, 1 Heisenberg models have a higher "hidden 
symmetry" compared with S — 3/2. Actually, there is 
another hidden symmetry for the spin 5 = 1/2 systems 
which is absent for spin-1. 



Xij Xij 



.-. At 



A^A., 



^N,N, 



-t/ 



where Ni = ClCi is the particle number on site i. Our 
analysis suggests that smaller spins have higher symme- 
try when the spin interaction is expressed in the fermionic 
representation. 

For 5=2, the Heisenberg Hamiltonian can be written 
as 



H=JY,^^ 



J 



S.^=-9ETr:(^JlV',-V'H) 



(».J) 



(so) 



',ti 



As in the case of lower spins, the terms Tr : {xpjltjjj ■ 
tfjjlipi) : and Tr : (iplipjtp.ipi) : have finite overlaps and 
are not completely independent of each other. For 5 > 2, 
it is not possible to represent the Hamiltonian in the 
above two terms alone. Quintet and higher multipolar 
hoping and pairing operators are needed to represent the 
Heisenberg Hamiltonian and we are not able to obtain 
any general expression. 

Summarizing, the Heisenberg model can be written in 
the fermionic representation as 

H = -^5][a(5)Tr:(Vl^,V']V.): 



(»j> 



+b{S)TY : (V'll^j • ^]lij,) : + ■ 



(11) 



where a(5), 6(5), • • • are parameters dependent on the 
spin magnitude 5 and • • • represents quintet and higher 
multiplet terms. In the following, we shall study in detail 
a Hamiltonian which keeps the first two terms only. The 



Hamiltonian can be considered as an effective Hamilto- 
nian for constructing trial ground state wavefunctions in 
a variational calculation. Notice that the effective Hamil- 
tonian is in fact "exact" up to 5 = 2 if a(S) and 5(5) are 
chosen properly. 



B. Path integral formalism and constraints 

In imaginary time path integral formalism, the parti- 
tion function is given by 

Z = Tr(e-^^) = f D^DxP^DXe'foL{i^.^'>^)'ir 



H 

(12) 



L = Y.\ [TrV^^(^- - ^\^oM\ ± ^(2^ - 1)A. 



where A^ is the Lagrange multiplier field introduced to 
impose the constraint Eq. ([4]). The +(— ) sign corre- 
sponds to the particle {N = 1) and hole [N — 25) 
representations. The Grassmann number t/i satisfy the 
antiperiodic boundary condition V'(O) = — '0(/3). 

Notice that the Lagrangian Eq. (|12l) is not invariant 
under the action of the internal symmetry group because 
of the non- invariant particle number constraint Eq. ^ , 
i.e. the internal symmetry group of the spin operator is 
not a symmetry group of the Lagrangian. In the follow- 
ing we will employ a trick to restore the symmetry of 
the Lagrangian. The integer- and half-odd-integer- spin 
models will be discussed separately because of the intrin- 
sic difference in their internal symmetry groups. 



1. Integer spins 

Our trick is to consider an average of all the partition 
functions with different particle-number constraint rep- 
resentations. In the integer spin case, the average is given 
by 



Y[ i [S{N, ^1) + S{N,- 25)' 



{i.r} 



= /i?^I?V^me-/ff(^.^Tr[^.a.^t]+^)rf, ^ 



TT i /'g*^[Tr(^,a,VJ)-2S+l] 
{i,r} 



(13) 



Notice that we do not make any approximations in de- 
riving Eq. (IT^ from Eq. p^ and the "new" partition 
function is a faithful representation of Heisenberg model 
except that it averages over all possible "particle" and 



"hole" representations of constraints locally. Similar 
ideas have been applied to generate the supersymnietric 
representation^ and the SU{2) representation of the t-J 
model-- . The fermion particle-hole symmetry is restored 
in this representation. The Lagrangian corresponding to 
Eq. (dSl) is 



L 



Ij: 



-H. 



TT[iP,{dr - jA,(7^)VI] - 2 In cos 



(25-l)A, 



This form of Lagrangian is invariant under the internal 
symmetry group U{l)<SiZ2 and the symmetry group be- 
comes almost a "gauge symmetry" of the new Lagrangian 
as we shall see in the following. We note that the La- 
grangian is complex because of the multiplier A. This 
problem can be solved by lifting the contour of integra- 
tion over A into complex plane via analytic continuation 
to find a saddle point in the imaginary axis^^. For conve- 
nience, we write Xi = iXi, and integrate Ai from —ioo to 
ioo (the saddle point of A^ is real) . Then the Lagrangian 
becomes 



i 

+H, 



TT[ipi{dT - Xiaz)ip. 



2 In cosh 



(2^-l)A, 



(14) 



Notice that A^ changes sign under an uniform Z2 gauge 
transformation, and Xij i ^ij exchange their roles under a 
"staggered" Z2 transformation where we have ax on one 
sublattice and I on the other sublattice. Notice also that 
the Lagrangian is not invariant under a time-dependent 

gauge transformation, because the term 2 In cosh - — ~ ' ' 



is not invariant. A similar situation occurs in the su- 
persymnietric representation of Heisenberg modeli-^. Be- 
cause of this restriction the internal symmetry group does 
not generate a complete "gauge symmetry" in the path 
integral formalism. 

Integrating out the fermion fields, we get an effec- 
tive action in terms of the Uij, V,y and A^ fields. The 
mean field values of the these fields are given by the 
saddle point of the effective action determined by the 
self-consistent equations U = (Uij), V = (V.y) with the 
mean field Hamiltonian 



Hm ^ JY. [a{S)[-{x*CJC, + A*CJC, + h.c.) 



M 



(^J> 



|A|2] + b{S)[-{v* ■ CJIC, + u* ■ C}lCj + h.c.) 



In cosh 



(2S'-1)A, 



(17) 



We shall now consider the path integral representa- 
tion in terms of the effective Hamiltonian Eq. (fTTj) keep- 
ing only the first two terms. The effective Hamiltonian 
can be decoupled by a Hubbard-Stratonovich matrix field 
(see Appendix B) 



U,. 



V,. 



i^li'j 



V.IV^, 




CJICo, and Vi 



CjlC,. Note that Ui,- and 



where u^ 

Wij form two sets of spin triplet operators, respectively. 

The Hamiltonian becomes 



H 



2 '^ 



{a{S)TT[ulu., 



(uyUj 



h.c. 



6(5)Tr[V^ • v., - (V., • VJIV-. + h.c.)]} {15) 



where x = {CJC)), A = (CJC,), u - (CJIC,), v 
(CjlCj), and A is determined by the condition 



(K 



2S+1 



25-1 , (2S'-1)A 
— - — tanh ^ — —. 



Notice that the averaged particle number satisfies 1 < 
(Ni) < 2S because we are mixing the particle and hole 
pictures in the constraint. If the ground state doesn't 
break the particle-hole symmetry, then A = and we get 
the half-fihing condition (N,) = {2S + l)/2. 



2. Half odd integer spins 

We first examine how the two particle number con- 
strains Tr(V'(Jz'(/'^) = ±(25*— 1) are transformed into each 
other by the SU{2) internal symmetry group. We con- 
sider the vector constraint Eq. ([5]) for S > 1/2 



The IJij,Vij and A^ fields define U{1)(^Z2 lattice gauge 
fields coupling to the fermionic particles which we shall 
call spinous in the following. Under the local gauge trans- 
formation ipi — >■ ipiWi, the temporal component and the 
spacial component of the gauge fields transform as 



X^<7, -^ W^{X,a,-—)W, 



dr' 



U,. 



1l 



-^ wju^.w,, 



IJ I 



(16) 



where 



M 



1 



25-1 



AU 



M, 



Tr(VaV'^) = (0,0,1)^, 



1 



25-1 
i 

25-1 



Tr(C'C^ + CC^), 
Tr(C'C^ - C&), 



M. - ^5^Tr(CCt - CCt). 



The Mx and My components indicate that the one- 
site pairings csc-s,- ■ ■ , cic_i are equal to 0, and the 
Mz component imposes the constraint that there is ex- 
actly one fermion per site. Notice that the first two 
components of the constraint is automatically satisfied 
if the third component is satisfied rigorously. Under the 
SU{2) transformation ip — > i/jW, the constraint becomes 
2^Tr(V'(7^t) ^ i?-i(o,0, 1)^, where WaaW^ = Rab<Jb, 
R is an SO(3) rotation. Since the constraint is invari- 
ant under an C/(l) gauge transformation generated by cr^, 
the continuum space formed by different constrains is the 
surface of a sphere SU{2)/U{1) = S"^. In particular, for 
W = 6*°'="'^/^, the north pole of the sphere which corre- 
sponds to the positive sign in Eq. ^ is transformed to 
the south pole, i?~^(0, 0, 1)^ = (0,0,-1)"^, correspond- 
ing to the negative sign in Eq. Q. 

Now we construct a Lagrangian which is invariant 
under the SU(2) internal symmetry group. In the 
path integral formalism, the particle number constraint 
2^Tr(i^,aV'l) = {Mt, Mf , Mf )^ = (0, 0, 1)^ can be re- 
alized by introducing a vector Langrange multiplier field 
A, 



S{MnS{M!)S{M:-l)^{^f 



dXfdXUXf 



exp{z[A^Afr + X^My + K{Mt - 1)]}. (18) 

The vector (0,0,1)^ changes into i^-^O^O, 1)^ = 
{rf ,71^ ^n'^Y' = ii under SU{2) gauge transformation, 
where n is an unit vector on the surface of the sphere 
S'^. As in the integer spin case, a SU{2) invariant path 
integral formalism can be obtained by averaging over all 
possible particle number constraints: 

= -^{^f fd^hJd^X,exp{i\i-{M,-n)} 

= i^f fd'X.,eMKK-M,-2ln^-^)} (19) 



then UijjSfij and \i form an SU{2) lattice gauge field 
coupling to the spinous. (Notice the sign difference in the 
matrix elements of U and V between integer and half- 
odd-integer spins.) Under gauge transformation tpi — >■ 
V'jWi {Wi £ SU{2)) the gauge fields transform as. 



A, -a ^ wJ{X,-a-—)W, 

and the Lagrangian Eq. (j20p is invariant only under time- 
independent gauge transformations as in the integer spin 



ill Xj 



case because of the non-invariant term (25' — 1) In ■ 

This is different from the S = 1/2 case, where the La- 
grangian is invariant under time-dependent gauge trans- 
formations. The mean field effective Hamiltonian corre- 
sponding to Eq. (f20| is 

i?™ = jJ2[^(S)[-{x*CJC,+A*C}C,+h.c.) + \x\^ 

+ \Af] + b{S)[-{v* ■ CJIC, + u* ■ CJIC, + h.c.) 

E25 — 1 ,r ,' , sinhA, 
^^(A-M,+ln^^) 



|uP + |vPl 



(21) 

where the parameters %, A,u,v are solved self- 
consistently by the mean field equations U — (Uij), 
V = (Vjj) and A is determined by 

(Mi) + -^(AcothA-l) =0. 
A^ 

The averaged particle number per site is again given by 
{Ni) = (25* + l)/2 if the ground state respects particle- 
hole symmetry (A = 0). 



where A^ = y/iXf)^ + (A^)2 + (Af )2. Introducing A, = 
iXi, we obtain a Lagrangian with "almost" SU{2) gauge 
symmetry: 



C. An important difference betw^een integer and 
half-odd-integer spins 



L ^ -Y^[Tr[ip^{dr -A,, •o-)V'|] - (25- l)ln^^^i — -] I* can be proved (see Appendix B) that for integer 



-H. 



(20) 



As in the case of integer spins, we replace H by the 
effective Hamiltonian Eq. dTTT) with only the first two 
terms. Introducing the Hubbard Stratonovich field as in 
the integer spin case (see Appendix B for details). 



^\\^, 






spins, Xjj = xlj, Aj, == -Ay, u^j = u^ and v^j = v^, 
whereas for half-odd-integer spins, Xji = xL? ^ji = ^iji 
Uji = — u,y and Vji — vt . Notice that the parity of the 
paring terms A^j and u^ are different for integer and 
half-odd-integer spins. These operators are central to 
the mean field theory, because the ground states are com- 
pletely determined by their expectation values. In partic- 
ular, for a mean field theory with A, x 7^ and u, v = 0, 
the mean field ground state is a BCS spin-singlet pairing 
state where the order parameter (Aij) has even parity 
for half-odd-integer spins, but has odd parity for inte- 
ger spins. We shall see how this important difference 
leads to different excitation spectrums between integer 



and half-odd-integer spin systems in the fermionic mean 
field theory. Similar result exists for states with A, x = 
and u, V ^ 0. In this case the mean field ground state 
is a BCS spin-triplet pairing state and the parity of the 
order parameters are reversed. 

We note that as in the S = 1/2 case, the mean field 
Hamiltonian should be viewed as a trial Hamiltonian 
for the ground state wavefunction of the spin systems 
after Gutzwiller projection. For the mean field theory 
with particle-hole symmetry the Gutzwiller projection is 
rather non-trivial. The state is a coherent superposition 
of states which allows sites with both one fermion and 
with 25 fermions. The numerical analysis of such a state 
is complicated and we shall not go into details in this 
paper. 



IV. FERMIONIC MEAN FIELD THEORY IN ID 
AND HALDANE CONJECTURE 



In this section we apply our mean field theory to the 
antiferromagnetic Heisenberg model in one dimension. 
We first consider the cases of spin 5 = 1 and 5 ~ 3/2 
where two versions of mean field theories based on dif- 
ferent methods of implementing the constraints will be 
discussed. The mean field results are summarized in sec- 
tion IV. C where the case of general spin 5 and Haldane 
conjecture will be discussed. 

To simplify our analysis we make use of the Wagner- 
Mermin theorem which asserts that a continuous sym- 
metry cannot spontaneously break in one dimension-^. 
Therefore we assume in our mean-field theory that the 
ground state is a spin singlet (spin liquid state) and the 
expectation values of (uy) and {vij) are zero, since the 
rotational symmetry will be broken otherwise. As a re- 
sult, we keep only the a(S) term in the trial Hamiltonian. 
For simplicity we shall also restrict ourselves to transla- 
tional invariant solutions of the mean-field theory in this 
paper. 



can be written as 

H,n = Xl^'=('^Ife^«'= + ^afe^'^'= + 4feCofc) 
fc 

^* 

- 2j[~^(Cs-fcCsfc - Ca-kCak - CQ-feCofe) + h.c] 
k 

+JN{\x\^ + I A|2) - N{^\ + In cosh ^), (23) 

where Xk = ^ — 2Jxcosk, A^ = 2iJAsinA: (since the 
phases of x ^nd A are unimportant in ID, we choose x 
and A to be real numbers) and N is the length of the 
chain. The mean field Hamiltonian can be diagonalized 
by a Bogoliubov transformation: 

Isk = UkCsk +Vkcl_f,, 



lak = UkCak - '^kc\_f,, 

7ofe = UkCak - «/cCo-fe- 



(24) 



The coefficients Uk,Vk satisfy the Bogoliubov-de Gennes 
(BdG) equations. 



EkUk = XkUk - A*kVk 
EkVk = -AfeUfe - XkVk- 

Solving the equations we obtain 



(25) 



Ek = V{\ - 2Jxcosfc)2 + (2JAsinfc)2, (26a) 

Ok Ok 

Uk = cos—, Wfe = zsin— (26b) 



where Ok is given by tan^^ — ^-^ and the diagonalized 
Hamiltonian is 

Hra = Y^ Ek{-fl,,-fsk + llklak + TofcTOfc) + Eq, 
k>0 

where Eo = EkiXk - Ek) + JN(x^ + A^) - 7V(|A + 

In cosh -1) is the ground state energy. Minimizing Eq, 
we obtain the self-consistent mean field equations at zero 
temperature. 



A. integer spin: S = 1 

We choose Htr to be the same as Eq. ([S]), i.e., a(l) = 
— 1, &(1) = 0. For convenience of discussion, we introduce 
two new fermions Cs , Cq which are the symmetric and 
antisymmetric combinations of ci and c_i. 



Csi = (cii +c_ii)/\/2 

Cat = {cii - C_ii)/\/2 



(22) 



The mean field Hamiltonian Eq. (IT7)) becomes completely 
decoupled in terms of Cs,Ca and cq. In Fourier space, it 






2N 



k 



Xk_ 
Ek 



. 25+1,-^ . Afe 
A = > sm/c , 



2N 



Ek 



'''^' N.)^'-^T. 



), 



Xk 



(27a) 
(27b) 



2N ^ Ek 

25-1 , (25-l)A ,_ , 
^tanh^ ^-(ZTc) 



with 5 = 1. The above equations are solved numerically 
where we obtain x = A = 3/4 and A = as shown in 
the first row of Table ID The averaged particle number 
of the spinous per site is 3/2, indicating that the ground 
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TABLE I: Solutions of two versions of mean fields for spin-1 
AF Heisenberg chain. E is the ground state energy per site 



and Egap is the excitation gap. 








(N,) X A A 


EiJ) 


^gap\^ ) 


(S. • Si) 


3/2 3/4 3/4 
1 0.7625 0.5671 0.7709 


-9/8 
-0.9032 


3/2 
0.3634 


3/2 
4/3 



state respects particle-hole symmetry. The ground state 
energy is Eq = — J(x^ + A^), and the excitation spectrum 
is flat (i.e. /c-independent) with a finite energy gap 3 J/2. 
The spin-spin correlation function is 



(S,; • S 



■i+r/ 



3 ^ 

p,q 


i(q-p)r 


Ep Eg 





Ep Eq 



and is nonzero only for r = and r = 1, with (S^ • S;) — 
3/2 and (Si-Si+i) = 3/4. The correlation function is zero 
for r > 2, indicating that the mean field theory describes 
a short-ranged Valence-Bond-Solid (VBS) state. 

We have also studied the mean field theory based on 
the Lagrangian Eq. p2]) for comparison. It gives rise 
to a dual of mean field theories with particle number 
constraints {Ni) = 1 and {Ni) = 2, respectively. There 
is a one-to-one correspondence between the solutions of 
the two mean field theories and we will only consider 
the "particle" representation {Ni) — 1. The mean field 
equations are the same as Eq. (P7)) except that Eq. (|27cp 
becomes 



{N^) 



2N 



a-E 



Xfc , 
Ek' 



1. 



The solution is also summarized in Table ID The excita- 
tion spectrum is gapped but no longer dispersionless. We 
note that the mean field solution Eq. (|27|) with particle- 
hole symmetry has a better ground state energy. 



B. 



half-odd-integer spin: S — 3/2 



For the spin-3/2 Heisenberg chain, we consider the ef- 
fective Hamiltonian Eq. (|2T|) with 6(|) = such that 
the grounds state is a singlet state. We shall also take 
a(|) — yI for reason we shall see later. In momentum 
space, the effective Hamiltonian becomes 



-Ei^li 



C_3_t.C3, 



fcCi^,) + /i.e.] + const., 



where Xfe = A 
fjA* CO! 
the following 



15 



Jxcosk and At 



iA„ 



^ JA* cos k. We shall set x and A to be real numbers in 



As in the 5 = 1 case, the Hamiltonian can be diag- 
onalized by Bogoliubov transformations and lead to the 
self consistent mean field equations at zero temperature, 

25* + 1 



X = 



2N 



.^cosfc(l-l^), 



25 + 1^ Afc 

A — > cosfc- 



2N 



(M+) 



Ek 

ft.: 

25* + 1 Y^ Afc _ (1 - A coth A) A+ 



N{2S~l)-^ r.k 



A2 



25+1 

N{2S~ 



E Xk 
., . E,, 



Xk _ (1 - AcothA)A, 



A2 



,(28) 



where S = 3/2, M+ = M^ + iMy, A+ = A^^ + i\y and 

Ek = ^{K - yJxcosfc)2 + |A+ - ^JAcosfcP (29) 

is the spinon dispersion. 

Solving the equations we find that A = (i.e. the 
ground state does not break particle-hole symmetry) and 
there exists infinite degenerate solutions for x a-nd A 
satisfying ypTJAp = 1.2732 (see Table II). This is 
a direct consequence of the SU{2) gauge symmetry for 
half-odd integer spins we mentioned in Section III.B.2, 
and all these solutions are equivalent. Notice that the 
ground state energy Eq = — ^^^(x^ + I^P) is the same 
as the expectation value of the Heisenberg Hamiltoian 
Eq. ([TU)) . which is why we should choose a(|) — y| in 
the effective Hamiltonian. As a result of the particle-hole 
symmetry (A = 0), the spinon energy dispersion Eq. (j29p 
is gapless at Fermi points k = ±7r/2. The spin-spin cor- 
relation is given by 



(S,: • S 



i+rj 



15 ^ 

47V2^' 

P:Q 


i{q-p)r 


Ep Eq_ 





Ep Eq 



and decays at large distance as (Si • Si+r) oc r~2 because 
of linearized spectrum around Fermi surface. This is also 
confirmed directly by numerical calculation. 

Similar to the S = \ case, we have also solved the mean 
field theory with particle number constraint {Ni) = 1 
or equivalently (M) — (0,0,1)"^. The solution is listed 
in Table [Til Since Az ^ in this case, the spin on dis- 
persion Ek = ^ \J{\zl2 — xcosfc)2 -I- |Acosfc|2 breaks 
particle-hole symmetry and has a finite gap over the 
whole Brillouin zone. The particle-hole symmetric solu- 
tion is also found to has a lower mean field ground state 
energy. 



C. Haldane's conjecture and Edge states 



Comparing the mean field energy dispersion Eq. (|26ap 
and Eq. (f^ . we find that the excitation spectrum of 



TABLE II: Solutions of two versions of mean fields for spin- 
3/2 Heisenberg chain. In tlie first row (x, A) means the real 



combinations satisfying ^y^ + Ap — 1.2732 




(iV>) X A A, E{3) E.aAJ) 


(S, -S.) 


2 X A -1.5198 
1 1.1490 0.7181 -1.2377 0.6732 


15/4 
45/16 



spin-1 Heisenberg model is gapped whereas the excitation 
spectrum for the spin-3/2 Heisenberg model is gapless in 
the mean field formulation with particle-hole symmetry. 
This difference between integer spin and half-odd-integer 
spin persists for any spin S*, if we consider the trial Hamil- 
tonian Eq. ^ or Eq. ^ with a{S) ^ 0, b{S) = 0, i.e. if 
we consider BCS spin-singlet ground state wavefunctions. 
The mean field equations for arbitrary S are the same as 
Eq. (I?7)) or Eq. ([25]) for arbitrary 5', and the dispersion is 
qualitatively the same as Eq. (|26ap or Eq. ((29|) . as long as 
we consider particle-hole symmetric mean field solutions 
which do not break translational invariance. The main 
effect of changing S in mean field theory is to change the 
number of fermionic spinon species. This result is con- 
sistent with the Haldane conjecture,™ which asserts that 
the integer spin Heisenberg chains have singlet ground 
state with finite excitation gaps and exponentially decay- 
ing spin-spin correlation functions, while half-odd-integer 
spin chains have singlet ground states with gapless exci- 
tation spectrums and power law decaying correlations. 
Our particle-hole symmetric mean field results agree well 
with Haldane conjecture. 

Haldane noticed that the time reversal operator sat- 
isfies T^ — 1 for integer spin and T^ = — 1 for half- 
odd-intcger spin, and this results in different Berry phase 
contributions from the topological excitations (skymion 
or instanton) in the path integral formulation given by 
gi27rSQ where 



Q = 



/dtdx — n ■ (dtu x dxTi) 
At: 



(30) 



is the Skyrmion number. For Q = odd integer, the 



Berry phase (e 



i2TrS\ 



is l'5 = 1 for integer spins and 



(— l)*? = — 1 for half-odd-integer spins. The difference in 
the ±1 factor between integer spin and half-odd-integer 
spin is the origin of Haldane's conjecture. The situa- 
tion is quite similar in our mean field theory. Notic- 
ing that Aij = CJCj is a singlet formed by two spin-S* 
particles. The Clebsch-Gordan coefficients implies that 

Aij for half-odd- 



iji 



-Aij for integer S and Aji 



integer S (see Appendix B). This sign or parity difference 
results in appearance of A sin k in Eq. (j26al) for S — 1 
and A cos k in Eq. ([^ for S — 3/2, which leads to differ- 
ent symmetries of ground state wavefunctions and differ- 
ent excitation spectrums between integer and half-odd- 
integer spin chains. 

The topological term Eq. ([50]) leads to exponentially 
localized edge states for open integer spin chains which is 
absent for half-odd-integer spin chains.^ This important 




FIG. 1: (Color online) The first Brillouin zone of the triangu- 
lar lattice. 6i = (27r, ^) and 62 ~ (0, ^) are two reciprocal 
vectors. The Brillouin zone is divided into three regular pieces 
by the three new bases u, v and w. 



difference between integer spin and half-odd-integer spin 
chains is also refiected in our mean field theory where zero 
energy Majorana edge fermions exist at the open bound- 
aries for integer spin chains which are absent for half- 
odd-integer spin chains. The existence of Majorana edge 
fermions for integer spin chains is a direct consequence 
of the odd-pairing symmetry for integer spin chainsi^ 
The details of the Majorana edge fermions is discussed 
in Appendix C. 



2D: 5 = 1 SPIN LIQUIDS ON TRIANGULAR 
LATTICE 



We show in the previous section that our mean field 
theory is able to capture the essential physics of spin- 
liquid states in ID antiferromagnetic quantum spin 
chains. In this section we shall apply our mean field 
theory to the 2D J1-J3 Heisenberg model on triangular 
lattice. We shall show that our mean field theory ad- 
mits new spin liquid solutions not explored before. Some 
results of this model have been reported in a previous 
paper^^, where particle number constraint is treated in 
the "particle representation" {{Ni) = 1). In this paper, 
we shall revisit this model with the particle-hole sym- 
metric constraint. The Hamiltonian of the J1-J3 model 
is 






Sj + Ji 



[id] 



s„ 



(31) 



where (i,j) denotes nearest neighbor (NN) and [i,j] 
the next next nearest neighbors (NNNN). For S = 1, 
there are two channels of decoupling the spin interaction, 
namely, S, ■ S, = -^ : Tr(V'JV',V]V'«) : or S, ■ S, ^ -^ : 
Tr{ipllipj-ip-Iipi) :. Since there is no Mermin-Wigner the- 
orem to protect us at zero temperature at 2D, we have 
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to keep both terms in the effective Hamiltonian Eq. ([TT]) . 
We shall choose b = 1 — a, with the weight a determined 
by minimizing the ground state energy. 



k=0 




J, 




y 


/ A 


i^" Ae 


a\ 


\ ~ 


A 


/ 




/4;r/3\ 


iSniy 



(a) 



(b) 



FIG. 2: (Color online) Two different pairing symmetries for 
the singlet pairing Aij. We indicate the relative phase of 
the nearest neighbor pairing (Aij) = Aie"^*^. (a), /-wave 
pairing. The dished lines are the zero lines of the pairing 
function (in lattice space or momentum space), (b), Px + ipy- 
wave pairing. There are no zero lines, the pairing function 
vanishes only at the F point. 



Next we introduce the mean field parameters xi, Ai, 
ui, Vi(for NN) and xs, A3, U3, V3 (for NNNN) to de- 
couple the Hamiltonian. We assume that u and v are 
parallel, so we only keep the z component of these vec- 
tors, and note them simply as ui^3 and vi^s (which are 
set as real numbers). As in 1-D case we consider here 
only spin liquid solutions that respect translation and 
rotational symmetries. Since the pairings have odd par- 
ity (Aij) = —{Aji), only p-wave and /-wave like pairings 
are allowed. Solving the mean field equations we find 
two kinds of solutions (see Fig. [2]), which have /-wave 
and Px + ipy-wetve pairing symmetries, respectively. We 
shall concentrate on these two kinds of solutions in the 
following. 

In momentum space (Fig. [ij, the mean field Hamilto- 
nian Eq. ([TT]) can be written as 



J^ 111, 



k 

k 

+ ^Ufc(4fcCu. - clifeC-ifc) 

k 

-^4(c_i_fcCife + /i.c.), (32) 



with 

Xfc = A - aZ(JiXi7fc + J3X3l2k), 

A* = iaZ(JiAiV'fc + J3A3V'2fc), 

Vk = -(1 -a)-Z'(Ji7fcWi + J372fcW3), 

ul = (1 - a)Z(Ji7fcWi + J372fcW3), 

where A is the Lagrange multiplier and 



i [cos kx 



Ik 



^k = i[sinA::r 



+ cos(-% + 4^)+cos( 



^/3ky 



-sin(--^ 



iff 



p+ip 
k 



[sin kr 






sin(^ + # 

2 ' )\- 



) + Bin(-^-^)], 

) 



The pairing symmetries can be verified easily by expand- 
ing the pairing terms at small k: A^ ex k^ik^ — 3fc^), and 
A^"" jpy ^ kx+iky. There are three lines of zeros for the 
/-wave pairing function as shown in Fig. [S] 

Introducing a vector Ck = (cife,cl-^_j,, coa:,Cq_j,)^, the 
Hamiltonian Eq. p2p can be written in a matrix form: 



Hrn — 2_^ Cf.HkCk 



J2^2^''~^''^' 



(33) 



where 



-u^. - Afc 



Hk^ 



I Xk+Vk 

-ul~^k -Xk + v^ 







V 



















Xk 


A. 


•> 


2 


K 


Xk 


2 


2 



Hk can be diagonalized by the Bogoliubov transforma- 
tion. 



Ak 



cos —Cik 



■ ^k ir,. t 

smye'-'^cLi^fe, 



-fe 2 

Dk == cos — co/c 



"Cik 



Ok t 
COSycLi_fc, 



Ok 
2 



sm—e'f"'cj 



O-fc' 



(34) 



where tan 0k == K±^ e*'"= = |"'-+^^ , tanOfc = -1^ 

'^ Xk ' |Mfc-l-Afc|' "^ Xk 

and e*'^'' = t^^ . The corresponding eigenvalues are given 
by Vk ± Elk and ±EQk, where 

Elk - V^x^ + JAfel^ + z.^, £;„, ^ V'x^ + IA.I^, 

and the self-consistent mean field equations are 

1 v^ r3 Xk Xk 1 



Xl,3 - 1^ 7 ,7fc,2fc[^ 
k 

Ai^3 = ^}^^M/c[^ 
k 



Xk 
Elk 



2E; 



Ok 



lAfcl 

2^;, 



Ofc 



1 



Wfe 



k 

^1,3 = 0, 

,.y 3, 1 , A 



(35) 
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where we have adopted the particle-hole symmetric con- 
straint in writing down the mean field equations. We 
have also checked the mean field solutions with constraint 
(Ni) = 1— and found that the excitation spectrums are 
qualitatively the same as those in the particle-hole sym- 
metric theory. The ground state energy is lowest for van- 
ishing wi,3,ui,3, i.e. a = 1, Uk = Vk = and the ground 
state is a spin-singlet. The excitations are characterized 
by three branches of fermionic spinous with Sz = 0,±1 
and identical dispersion Ek — \/Xk + l^fcP- 
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FIG. 3: (Color online) Phase diagram for the J1-J3 Heisenberg 
model with Ji + J3 = 1. We use the constraint {Ni — |) = 

^tanh-|. The red dot line indicate energy per site for the 
Px + ipy state and the black square line is for the f-wave state. 
Dirac nodes exist for the f-wave state. The number of nodes 
is 6 when Ji is dominating and 24 when J3 is dominating. 

Our mean field theory for the J1-J3 model contains 
two regimes of spin liquid states for both / and Px + 
ipy pairing symmetries as a function of J1/J3. A first 
order phase transition occurs between the two regimes at 
J1/J3 "^ 1 (see Fig. ^. When Ji dominates, the spin 
liquid state is characterized by xi,3 7^ a-nd Ai^3 7^ 
(consequently (S^ • S^+i) < and (Si • 8^+2) < 0); while 
when J3 dominates, xi = ^1 = and X3 7^ 0, A3 7^ 
(consequently (Sj • S,+i) = 0, (S^ • 8^+2) < 0). The 
Px + iPy states remain lower in energy in both regimes. 

The /-wave pairing solution respects particle-hole sym- 
metry and has (Ni) = |. The excitation is gapless with 
several Dirac cones in the first Brillouin zone. A cut 
of the spinon dispersion is shown in Fig. 2] where the 
particle- hole symmetry is obvious. The mean field solu- 
tion with {Ni) = 1 has similar properties, except that 
the particle-hole symmetry is lost and the position of the 
Dirac nodes are shifted.— 

For the Px + ipy-wave pairing, we find that the ground 
state breaks particle-hole symmetry and A 7^ 0. The ex- 
citation spectrum is fully gapped. Similar to integer spin 
chains, we find that the Px + ipy-wave state is topologi- 
cally nontrivial (see also Appendix C). The solution has 
|A| < Z\JiXi + J3X3I, so that Xfc can be either positive 
or negative, depending on k. Therefore the Bogoliubov 

' jAfc ReAk Xk \T 



(which is S"^) covers the Brillouin zone (k space, which 
is also S"^) at least once. In other words, the topological 
(Skyrmion) number of mapping from the k space to the 
spinor space is nonzero m ^ Q. However in the vacuum 
where the spinon density is zero, A is very big and Xk can 
only take positive values, which gives a zero topological 
number;^ Since the bulk and vacuum belong to different 
topological sectors, the boundary defines a domain wall 
between the two phases, and there should exists gapless 
(chiral) Majorana edge states in the Px + iPy state follow- 
ing the analysis of Read and Green. ^^ In this sense, the 
Px + ipy-state describes a time-reversal symmetry break- 
ing topological spin liquid. 




FIG. 4: (Color online) The spinon dispersion the f-wave pair- 
ing state with constraint {Ni — |) = i tanh -^ for Ji — 0.6, 
J3 = 0.4. The figure on the right shows the dispersion along 
the line linking two K points of the Brillouin zone (see the 
blue lines on the left). The negative energy part is filled and 
the excitation is gapless with 6 Dirac nodes in the first Bril- 
louin zone. The dispersion for the mean field solution with 
{Ni) = 1 is similar to the above one except that the locations 
of the nodes are shifted a little bit. 



spinor space described by the vector (- 



By 



Ek 



Next we consider the specific heat and spin suscepti- 
bility for the /- and Px + iPy states. For the Px+ipy-y^&'ve 
pairing states, since the excitations are fully gapped as 
in BCS superconductors, the specific heat and spin sus- 
ceptibility will show exponential behavior^^ at low tem- 
perature. For the /-wave pairing states, both the specific 
heat and the spin susceptibility show power law behavior 
due to the Dirac cone structure (See Fig. 01). The energy 
for the system '^^ E — i^j. ^^1'^'' in mean field theory. 
At low temperature, the specific heat is dominate by the 
low energy excitations near the nodes of the fermi surface, 
which is given approximately hy Ek = vpk^ where vp is 
the fermi velocity (the anisotropy of the dispersion can 
be removed by a re-scaling of momentum and will not af- 
fect the result qualitatively). The magnetic specific heat 
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is given by 

dE d f'^ vpk 

^M — T^^ — -J-;^^ / 1 2'Kkdk 

dT dTj„ e^+1 



2TTx'^^rdx 



-3A 

ex uttT /vp (36) 

where n is the number of Dirac cones. The magnetic 
susceptibihty can be calculated from the linear response 
theory with 



Xz{kM = / dV /drr,(5,(r,T)5,(0,0))e*('^^^'='-) 



2 



/" Sr I dre*("^-*='') ^ e*(«'- 



q)r 



x[T,(cI(g,T)ci(g',T)cl(p,0)ci(p',0)) 
-r,(cl((z,r)ci(g',r)c^_i(p,0)c_i(p',0))] 

(37) 

The static susceptibihty is given by Xz ~ 
hmfe_j.oX2(fc,0) oc -^. The detailed calculation is 

given in Appendix D. 

Recently, a spin-1 material NiGa2S4 was discovered, 
where the S ^ \ Ni^^ ions form a triangular lattice 
with antiferromagnetic (AFM) interaction. The system 
has a Curie- Weiss temperature ~ 80K, and is found to 
exhibit no conventional magnetic order down to 0.35K. 
Several plausible ground states have been proposed for 
this system^i^r— including an antiferro-nematic (AFN) 
staters, a ferro-nematic state,— and a KT phase driven 
by vorticesi^ Spin liquid is also proposed to be a plau- 
sible ground state. In our mean field theory, the /-wave 
pairing spin liquid has coherent gapless excitation with 
specific heat scaling as T^ at low temperature which is 
consistent with the experimental result.— Furthermore, 
the Knight shift data in Ref. [l^ shows that the intrin- 
sic susceptibility decreases with decreasing temperature 
(but the very low temperature data are missing.—) This 
matches qualitatively with our mean field result Eq. ([37| . 
We therefore proposed that the /-wave pairing spin liq- 
uid state is a plausible ground state for the material. 



VI. COMMENTS AND CONCLUSION 

Sumnrarizing, we show in this paper how the fermionic 
representation for spin S* = 1/2 systems can be general- 
ized to spin systems with S > 1/2. The symmetry group 
of the spin operator is SU{2) group for integer spin and 
U{1)®Z2 group for half-odd-integer spin. Different path 
integral formulations and mean field theories are devel- 
oped corresponding to different ways of handling the con- 
straints. In ID, we show that the particle-hole symmetric 
mean field theory for S — \ and S — 3/2 spin chains are 



consistent with Haldane's conjecture, and we argue that 
the difference reflects a fundamental difference between 
integer and half-odd-integer spin chains. We also study 
2D spin-1 AFM on triangular lattice where we find two 
spin liquid states, one is a gapless /-wave spin liquid and 
the other is a topological Px + ipy spin liquid state. We 
propose that the gapless /-wave spin liquid is a plausible 
ground state for the material NiGa2S4. Our approach 
can be applied to any other spin models and provides 
a new approach to spin liquid states for S > 1/2 spin 
systems. 

We note that because of limitation in scope we have ad- 
dressed only a very limited number of issues in the study 
of spin-liquid states in this paper. Within the mean- 
field theory we have restricted ourselves to spin liquids 
solutions of Hamiltonian pip with spin rotational sym- 
metry (ID) and lattice translational symmetry (ID and 
2D). We note that solutions which break spin rotational 
symmetry exist in our theory which may serve as ground 
states of Hamiltonian PTJ! ^'^. For example, we show in 
Ref. [2^ that state with long-ranged magnetic order exists 
as lowest energy state of the mean-field theory of the Ji- 
J3 model and anti-ferro nematic order may exist when 
the bi-quadratic {K{Si ■ Sj)^) spin-spin interaction exists 
in the Hamiltonian. For simplicity we only considered 
the mean field states without breaking the translational 
symmetry. We note that states that break translational 
symmetry (such as dimerized states) are believed to be 
ground states of some ID or 2D spin models. ^^ 

Another very important issue is whether the spin liq- 
uid states we find are stable against gauge fluctuations. 
In 1-D, the gauge fluctuations can be removed by a time- 
dependent Read-Newns gauge transformation and thus 
have no effect to the low energy properties^. The situa- 
tion is very different in two or higher dimensions where 
the stability of the mean-field state depends on dimen- 
sionality and the (gauge) structure of the gauge field fluc- 
tuations. For instance, Z2 spin liquid is believed to be 
stable at 2D2£ and gapless U{1) Dirac fermionic spin liq- 
uids is stable at 2D in the large- A^ limi1^. In our case of 
spin liquid solutions for the S = 1 Heisenberg model, the 
C/(l) gauge fluctuation is gapped by the spinon-pairing 
term via the Anderson-Higgs mechanism. The ground 
state is stable against the low energy gauge fluctuations 
which is described by an effective Z2 gauge theory. We 
note that the emergence of gauge-field structure is a di- 
rect consequences of particle number constraints and a 
reliable answer to the question of stability of spin liquid 
states can be obtained only if we can handle the parti- 
cle number constraint reliably. Gauge field theory can 
only handle long distance, low energy gauge fluctuations 
and a more satisfactory answer to the question of sta- 
bility of the spin liquid states can be obtained only af- 
ter the Gutzwiller Projection wavcfunctions arc studied 
carefully. 
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Appendix B: proof of some algebraic relations 

First we introduce a (25* + I) x (25* + 1) matrix B, 
which is essential in our proof: 



B = 



1\ 



(Bl) 



Appendix A: Symmetry group of the spin operators 
and the pairing parity 

We shall look for allowed transformations V' — ^ ipW, 
which keep the spin operator Eq. ([3|) invariant. At first 
glance, this condition is satisfied as long as W G U{2). 
However, it is not true because Cm's and cj„'s are 
not independent. By checking each row of rp directly, 
it is easy to sec that (cm, (— l)'^^™cL„j) transforms 
to (c,„, (-l)^-"cl„jW^, and (4, (-l)^-™c_™) ^ 
(cJ„,(-l)'^~"c_™)VK*. Noticing that 



(cL,(-i) 

obtain 



S-m, 



= (-i)^-"(c-m,(-i)^-'"4K, 



or 



(-1)^-™(c_„„(-1)^-'"cLK-^ 

(c_„,(-l)^-'"4)^ 
(c_„,(-l)^-™r4)a.l^V,. (Ala) 

On the other hand, (c_m, (— l)'^^™cj„) trans- 
forms to {c-„,,{~lf+'^cl^)W. Noticing that 

(c_,„,(-l)^+-ct„) = (c_„.,(-l)-^-™4J = 
(c_„,(-l)^-™4)(a,)-2S and (a,)-2S = (a,)^^ 

we obtain 

(c„™,(-l)^-'"ctj(a.)2^^ 
(c_™,(-l)^-'"cL)(a.)^^W^. 
or 

(c_,„,(-i)5-™4j^ 

(c_™, (-l)^-'"4)(a,)2SW^(a.)2S. (Alb) 
Eq. (jAlaP and Eq. (jAlbl) impose the condition 

(TxM^V, = (a,)2^W^(a,)25. (A2) 

For integer spin, Eq. (IA2p becomes W* = a^Wa-j., 



which results in W 



or W 



'a^, i.e. W e 



U{1)®Z2, where U{1) = {e"^' } is the usual gauge sym- 
metry group and Z2 = {\ax\ is the particle-hole sym- 
metry group. Notice that the two groups don't com- 
mute, thus the total symmetry group is a semidirect 
product of the f/(l) group and the Z2 group U{1)®Z2 = 
{e"'''^,a^e"''^ = e-'^-^a^-.e e R}. For half-odd-integer 
spin, Eq. (|A2p implies W* = ayWay, which is equiva- 
lent to detVF = 1 since W^W = 1. This indicate that 
W G SU{2). Thus the symmetry group is SU{2) group. 



The matrix is the CG coefficients combining two spin-S to 
a spin-singlet state, i.e. |0, 0) = Bjnmi\S^m)i\S,m')2 — 
C^C\vac). Thus C can be expressed as C = BC* — 
_B(c^,c^_-^, • • • T<y_gY' ■ Notice that for integer spin 
B~^ — B while for half-odd-integer spin B^^ — —B. 

The following are some common properties of the B 
matrix. Suppose R is an rotation operation, and D{R) is 
its (25* -1-1) X (25* -1-1) dimensional matrix representation, 
then we have RC = D{R)C, and 

RC = D{R)C = D{R)BC*. 

On the other hand, RC* = D{R)*C*, so 

RC = R{BC*) = BRC* = BD{R)*C*. 

Comparing these two, we have D{R)B = BD{R)* , or 

B-^D{R)B^D{R)*. (B2) 

This means that B is the transformation matrix which 
transfers the representation matrix of a rotational oper- 
ator to its complex conjugate. In general, D{R) can be 
written as D{R) — e*'-'^, where R are the SU{2) group 
parameters (real numbers) for the operator R, and con- 
sequently D{R)* = e"*^*'^. From Eq. ^^ we have 

B-'D{R)B - B-ie'i-^S - e''^"i^-^ 

-jI* R 



D{Ry 



So we have 



R-^IB = I* = I^. 



Noticing B ^ = ±B, above equation is equivalent to 

Now let us consider the spin algebra. The spin inter- 
action is decomposed by the following two operators 



^Hj 



C]C, C]C 



c]c 



c]c. 



V-Ii^, - \A 



3 ^i ^1 



C\\Cj C\\Cj 
CjIC, CjlC, 



(B3) 
(B4) 



We define Xrj = CJCj, A,, = CjCj and u.^ = CJICj, 
Vij = CJICj. The following discussions distinguish 
integer- and half-odd-integer- spins and prove the results 
in section III.C. 
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integer spins {B^ = B^ = B^^ = B) 
The second term on the first row of Eq. (|B3p can be 
written as 



4/ 



-i*\T 



CJC, = CJBC; = ~{C}BC* 

= -c]bc; = -c]c, = -aI 

In second step, we have used the fact that the transverse 
of a scalar operator is equivalent to reversing the order of 
the constituting components. Commuting two fermionic 
operators we obtain a minus sign. On the other hand, 



^tz-.U 



-t, 



qc,, soA,, 



-A, 



The second term on the second row of Eq. (|B3|) is 



cjCj = cf BBC* = -{ct qy = -qa = -x. 



,t 



It is obvious that 



J 



To simplify Eq. (|B4I) . we will use the following rela- 
tions: BIB = -I^, Bl = -l^B and IB = -Bl^ . The 
second term on the first row of Eq. (jB4[) is 

cliCj = cliBq = -qBi^q = {qBi^qf 

= C]lBC: = C]lC. = {CllC,)^ ^ nl (B5) 

On the other hand, above operator can be written into 
another form, CJICj — {qiCi^ — u'-^. Comparing with 



Eq. (jB5|) we get u^ = Uji. 

The second term on the second row of Eq. (|B4[) is 



CJICi = CjBlBC* 



-cTi^c* 



= {cji^qf ^c]ic,^w. 



(B6) 



And the relation vj^ — Vji manifests itself. 

half-odd-integer spins {B'' = B^ = B-^ = -B) 
Repeating the above procedures, it is straightforward 

to show that 



etc, 
chq 



-ut cllQ=vl 



with 



Xji — Xijy Aji — iXij^Uji = — Uy, Vji 



(B7) 



Appendix C: Edge states in open Integer spin chains 

It is known that integer spin antiferromagnetic Heisen- 
berg model has free edge states with spin magnitude- 
S/2.^- Majorana fermion edge states also exist in our 
mean field theory under open boundary condition. Our 
discussion will follow the argument of Ref. [23 for 2D 
pairing fermions. Firstly, we point out the existence of 
two topologically distinct phases under periodic bound- 
ary condition and then we discuss the zero energy edge 
state solution at an open boundary. 

In our mean field theory, all the properties of the 
ground state are completely determined by Uk and 



Vk (or Xk and Afc). Noticing that Uk and Vk obey 
Iwfcp -I- |wfcp = 1, so they can be viewed as a spinor. 
This spionor is equivalent to a vector (u* ,v*)a{u,v)^ = 

{0, sin 9k, cos 9 k)'^ = {0,-[^:^)^ where we have used 
Eq. (j26bp . Obviously, the spinor {u, v) spans a S^ space. 
Recalling that the first Brillouin zone is also 5^, the 
functions Uk,Vk describe a mapping from S'^(k space) 
to S"^ (spinor space). The mapping degree is character- 
ized by the first homotopy group 7ri(S'^) = Z. When 
—2Jx < A < 2Jx (which is the case for our mean field 
theory) , the map is nontrivial because every point of the 
spinor S^ is covered once, in other words, the mapping 
degree is to = 1. This topological number defines a phase, 
we call it A phase. When |A| > 2Jx, the mapping is topo- 
logically trivial since Xk is always positive(or negative) 
and the lower half circle(or upper circle) of the spinor S^ 
is never covered, so the mapping degree is to = 0. We 
call this region B phase. Since the mapping cannot be 
smoothly deformed from 7TT, = 0torn=l,a topolog- 
ical phase transition occurs at |A| = 2Jx- We will see 
that the existence of zero energy edge states is tied with 
a phase boundary between the bulk (A phase) and the 
vacuum (B phase). 

Next we consider an infinite chain with a single edge. 
We assume that the edge is located at a; = 0, and x < 
is the vacuum where the spinon density is zero. The 
boundary can be described by a potential V{x) that is 
large and positive at x < so that the wavefunction of 
the spinon vanishes exponentially at a; < 0. Equivalently, 
we can assume a position dependent lagrange multiplier 
term A(a;) which becomes very large and positive(so that 
the condition |A| > 2Jx is satisfied) at a: < 0. This 
implies that the vacuum belongs to the B phase and the 
boundary is a domain wall between A phase and B phase. 

We now study the edge states in the continuum approx- 
imation. It is convenient to introduce iJ,{x) — 2Jx— A(a;) 
such that fi is positive in the bulk, becomes zero near the 
edge and turns negative outside the edge. Expanding Xk 
and Afe to first order in k near k — 0, we obtain Xk ^ — M 
and Afc ^ 2iJAk. To see the effects of the domain wall, 
we consider the BdG equation Eq. (P5|) . which can be 
written in position space as 



(CI) 



du 

at 


^ ^ A dv 

= -nu-2JA—- 

ox 


dv 


r. ^ ^ du 
= iiv + 2JA — , 
ox 



where we have replaced E and fc by i^ and 



' dt 



"dx^ 



respec- 



tively. The BdG equation is compatible with u{x, t) = 
v{x, t)*, the spinor (u, v) satisfying this relation describes 
Majorana fermions. There is a normalizeable zero energy 
solution for the above equations. Putting u — v, then 



2JA— = 
ox 



-flU, 



and the solution is 



u{x) — v{x) —(X e 2j/ 



■ f^ ^{x)dx 



15 



Notice that the BdG equation Eq. ((CT 



Eu 

Ev 



/if 



^lu - 2 JA 
2JA 



dv 

dx 
du 



dx 



at -E ^ admits no normahzable bound solutions. This 
shows that the zero mode is the only plausible bound 
state solution. 

In the case of a pair of boundaries at x = and x = W, 
we can replace the two first order equations by a second 
order equation, 



d^ 



d^ 



E\u ±v)^ [^i' - (2JA)^-^ T (2 JA)-^](u ± v) 



The above discussion is applicable to our mean field 
theory for 5 = 1 open antifcrromagnetic chain. The dis- 
cussion can be generalized to larger integer spins 5* > 1 
if we adopt the effective Hamiltonian Eq. (|23|) . 

Unlike the 5=1 case, our S = 3/2 mean field theory 
produces a mapping from fc-space S^ to the (u, v) space 
S^ which is topologically trivial (essentially because of 
the even parity of the pairing term) , and the Majorana 
edge states don't exist anymore. This is in agreement 
with the result that no exponentially localized edge state 
exists for open half-odd-integer spin chains. The ex- 
istence of power-law localized edge states in half-odd- 
integer spin chains^^ is more subtle and is not produced 
by our mean field theory. 



dx'^ 



dx 



and the energy E is not zero for finite W . However two 
bound state solutions with energy going to zero expo- 
nentially with increasing W exist. The solution u -I- w is 
centered at a; = and u — v\s centered dX x = W . Notice 
that there is only one solution for the above BdG equa- 
tions, and the two modes m ± w are Majorana fermions. 



Appendix D: Calculation of the susceptibility 

Substituting the Bogoliubov transformation Eq. ([M)) 
into the Kubo formula Eq. ([37)) . we obtain 



J 



-CqSq[GA{q, ~r) + GB{-q,r)] x C^+kSq+k[GA{q + fc, r) + Gei-q - k, -r)]] 
= ^ E [K^Aiq, in) - S^Gsi-q, -in)][-Cl^^GA{q + k,tn + lio) + S^^^Gsi-q - k, -tQ - lio)] 

^ q,in 

~CqSqCq+kSq+k[GA{q,in) + Gsi^q, -in)][GAiq + k,iQ + iuj) + Gsi-q - k, -iQ - iuj)]] 

a a 

with the notations Cq — cos -j and Sq = sin -j. Using the results GA{k, iuj) = Gsik, ioj) = G{k, iuj) and G{—k, iuj) 
G(k,iLu) — ^^_^ , we get the static susceptibility 



Xz = limxz{k,0) 



fc-s-O 



lim -- Y, [-{Ccfiq+k + SqSq+kfG{q, in)G{q + k, in) + {CqSq+k - SqCq + kf G{q, in)G{q + fc, -l^!) 



fc->o PV 



q.iO, 



= hm-V 



fe-i-o V 






£q+k — £q 



= -2 / d^q 



,2driF{eq) 



de„ 



(Dl) 



which is the Pauli susceptibility for free fermions. At low temperature, the excitation spectrum can be approximated 
by Dirac cones, so we have 



-2 / 2^q 

/o VpOq 



dn.isq)^^ 



— / nF{£q)dq ex -^ 
vf Jo vj. 



(D2) 
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